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Abstract 

We investigate the behavior of the entanglement entropy of a confining gauge theory near 
cosmological singularities using gauge/gravity duality. As expected, the coefficients of the 
UV divergent terms are given by simple geometric properties of the entangling surface in the 
time-dependent background. The finite (universal) part of the entanglement entropy either 
grows without bound or remains bounded depending on the nature of the singularity and 
entangling region. We also discuss a confinement /deconfinement phase transition as signaled 
by the entanglement entropy. 
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1 Introduction 

There has been considerable interest recently in the entanglement entropy of quantum field theories. 
Entanglement entropy is a measure of long range correlations in the system which has proven useful 
in a variety of applications including identifying exotic ground states [U |2] and signaling phase 
transitions. In particular, it has been shown O HJ |5] that a confinement/deconfinement phase 
transition, which is usually studied at finite temperature, can be detected in the zero temperature 
ground state by studying the entanglement entropy for regions of various sizes. 

We explore another potential application of the entanglement entropy: probing cosmological 
singularities. Understanding past spacelike ("big bang") singularities is a long-standing goal of 
quantum gravity. While much progress has been made in understanding static timelike singularities 
in string theory, spacelike singularities remain mysterious. There have been various attempts to 
resolve spacelike singularities using gauge/gravity duality (see, e.g., (6j [7J [8j [9]) but we will be less 
ambitious here. Rather than attempt to resolve cosmological singularities, we will endeavor to 
understand one of their effects on a quantum field theory. We compute the entanglement entropy 
for a confining gauge theory in a class of fixed (classical) cosmological backgrounds. As we will 
discuss, the choice of a confining gauge theory has technical advantages for the calculation, but it 
is also of particular physical interest since it models an early-universe quark-gluon plasma. The 
quark-gluon plasma is expected to feature a mass gap and a confinement/deconfinement phase 
transition at around 150 MeV [10J. 

Given a state of a quantum field theory and a region of space, 7Z, the entanglement entropy 
measures the quantum entanglement of the part of the state in 7Z with the rest of the state in TZ, 
the complement of TZ. It is defined in terms of the von Neumann entropy of the reduced density 
matrix on 7Z: 



where p^ = tr^p. While See is a divergent quantity due to the short-wavelength modes across 
dTZ and must be regulated, it is generally possible to extract its universal (UV cutoff-independent) 
component. 

Calculating See in strongly coupled theories using the field theory formalism, however, is 
notoriously difficult. We will therefore calculate See using the more tractable holographic approach 



See = -Tr (p^lnp^) 
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instead. The original Ryu-Takayanagi formulation [11 J of holographic entanglement entropy states 
that the entanglement entropy of a region 1Z in a field theory dual to some static asymptotically 
locally AdSd+i bulk spacetime is proportional to the area of a minimal surface homologous to the 
entangling region on the boundary: 

Area(7^) 

SEE ~ 7^ (L2) 

where 77^. is the minimal area surface in the bulk and G^ 1 ^ is the [d + l)-dimensional Newton's 
constant^ In this proposal, the classical bulk spacetime is static, so there is a canonical time 
slicing; the region 1Z and the corresponding minimal surface 7^ both lie on the same constant 
time slice. The more general, covariant Hubeny-Rangamani-Takayanagi prescription calls for the 
minimal surface 7^ to be replaced by an extremal surface; if several exist, the one with minimal 
area is chosen |12| . Since we work in cosmological backgrounds, we will employ the covariant 
proposal here. 

There have been earlier studies of entanglement entropy in time- dependent backgrounds, e.g., 
[1 ^ [H I IIB ] fI5 | 117 1 [TS j 11^ BU j [2T| but to our knowledge, none have explored the behavior near 
singularities^] 

The area of any bulk surface which extends all the way to the boundary at infinity diverges, 
and this divergence obeys [24J: 

A / \ °d-2 . Q-4 . / Co log (7) H d is even , . 

Area(7 R ) = —j-r + —7— H h< n_3j 

v ' e d ~ 2 e d ~ A \ c H a is odd v ; 

where L is the length characterizing TZ, e is a UV cutoff, and the dots in the brackets represent 
terms that vanish as we take e to zero. If a boundary-covariant cutoff is used, the coefficients q of 
the divergent terms are integrals of geometric quantities. In particular, the leading coefficient is 
proportional to the area of the boundary of the entangling region, &R. We will work with d = 5 
and determine the C\ coefficient as well as the UV cutoff-independent Cq as functions of time in a 
class of cosmological backgrounds. 

We treat the bulk geometry classically, but expect stringy and quantum corrections to be 
important very close to the singularity. We will stay away from this poorly understood region and 
only consider the entanglement entropy for times where the extremal surface stays in a region of 
the bulk where classical general relativity is valid. 



2 Cosmological backgrounds and extremal surfaces 

We wish to calculate the entanglement entropy for a field theory in a Kasner3 + i x S 1 background: 

ds 2 hdy = -dt 2 + t 2pi dx 2 + t 2p2 dyl + t 2p3 dy 2 2 + d6 2 (2.1) 

^^Note that this proposal admits quantum and stringy corrections when the bulk spacetime is not taken to obey 
classical general relativity. 

2 Thc interesting recent paper |23| considered the contribution to entanglement from the region inside black holes, 
but the extremal surfaces were bounded away from the singularity in the bulk. 
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where three of the spatial directions are periodically identified: 



9 ~ 9 + L 
Vx ~ Vx + Li 
Vi ~ Vi + ^2 



and the exponents £>j obey the Kasner conditions: 



5>? 



(2.2) 



These last conditions ensure that the metric is Ricci flat. This spacetime describes a homogeneous 
but anisotropic cosmology and generically has a curvature singularity at t = 0. The constant t 
surfaces are flat, but space expands at different rates in different directions. Since there are two 
conditions on the three exponents pi, there is a one parameter family of solutions labeled, e.g., 
by Pi which can range from —1/3 to 1. Generically at least one Pi is negative and the rest are 
positive, so one direction contracts while others expand. The only exception is when two of the 
exponents vanish and one recovers the Milne solution, which is locally flat. This is the only case 
in which t = is not a curvature singularity. 

We choose our entangling region, denoted TZ(tb, xq), to be a strip which is extended in all spatial 
directions except one, on a fixed time slice t = t 6 : 



n(t b ,x ) 



x G [-Xo,X ] 

9e[0,L ) 

yi e [o,LO 

2/2 e [0,L 2 ) 

t = t b 



(2.3) 



It is the goal of this paper to calculate the entanglement entropy of 7£(i&, xq) holographically as a 
function of t b . To do so, we must first find a bulk dual to the Kasner 3+1 x S 1 background. 

We start with the AdSg soliton metric [25 l l26 | . obtained by analytically continuing the timelike 
direction and one spatial direction from the planar AdS black hole: 



u * soliton 



(1 - z 6 ) d9 2 + j] liV dx ii dx v + - 



dz 2 



(2.4) 



where the AdS radius and black hole horizon radius have both been set to one. The 9 circle 



smoothly caps off at z — 1, provided we set L 



An 
5 ' 



i.e., we periodically identify 9 



In 

5 



+ 9. 



This solution describes a confining vacuum in the dual field theory which lives on IR 3 ' 1 x S 1 . One 



can change the boundary metric by using the fact that replacing r/^ in (2.4) with any (3,1)- 
dimensional, Ricci-flat solution of Einstein's equations still results in a bulk metric which is an 
exact solution of Einstein's equations with a negative cosmological constant. Replacing rj^ with 
Kasner yields the (5 + l)-dimensional Kasner-AdS soliton: 



r/s 2 
a6 KAS 



j2 (r^5 - df2 + t2pidx2 + t2p2d y 2 + t2pzd yl + ( x - z ") de2 ) 



(2.5) 
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This solution has a spacelike singularity at t = which extends across the entire spacetime. 

We could have started with AdS in Poincare coordinates, and replaced the Minkowski metric on 
each constant radial surface with Kasner. The result is again a solution of Einstein's equation with 
negative cosmological constant. This would be dual to a nonconfining gauge theory on the Kasner 
spacetime. However, in addition to the singularity at t — 0, this solution has singularities on the 
Poincare horizon which appear to be unphysical. These unphysical singularities are removed in 



(2.5) since the radial direction is capped off and there is no Poincare horizon. Another advantage 



of starting with the AdS soliton is that since it is a confining geometry, it can be used to probe 
the confinement/deconfinement phase transition. 

The proposal in jT2] dictates that entanglement entropy be calculated from the area of an 
extremal surface homologous to the entangling region. We use the symmetries of the surface and 
the bulk metric to parametrize the surface in terms of the distinguished boundary coordinate x: 

X M = (t, x, 6, 3/i, jte, z) = (T(x),x, 6, y 1: y 2 , Z{x)) (2.6) 

The induced metric on the surface 7 a fc = d a X' J 'dt ) X u g lll ,, where g^ u is the bulk metric, gives rise to 



the area functional for (2.6): 



A 



A7rLiL 2 



dxT ^iw~ [(1 " z[xf) {T{xfP1 " T ' {x)2) + z ' {x)2] 



(2.7) 



The set {T (x) , Z (x)} which solves the equations of motion generated by (2.7) and obeys the 
following boundary conditions 

Z (x = ±xq) = 
T(x = ±x ) = t b 
describes the extremal surface that enters into ( 1.2[ ). 



(2.8) 



We pause here to note that the integral in (2.7) diverges and must be regulated. It is conven- 
tional to impose a UV cutoff some small distance away from the boundary z = e, which is the 



same UV cutoff that appears in the expansion (1.3). The integral in (2.7), however, is over the 



variable x, so we define the UV cutoff to be the value z = e such that Z(xq — 5) = e. We then 
integrate up to x = Xq — 5 instead of up to Xq. 



2.1 Finding the Extremal Surfaces 



The equations of motion generated by (2.7) are a set of two coupled second order ordinary dif- 
ferential equations for which an analytical solution is not known. The general solution depends 
on four parameters, but the local boundary conditions (2.8) reduce this to two. We will generate 
a 2-parameter family of approximate solutions to the equations of motion by truncating a set of 
series expansions, and then use numerics to determine those two parameters for different values 
of the power pi and the time slice of the entangling region tb- We now describe this procedure in 
more detail. Since it is valid in arbitrary dimension, in this subsection we work in general d. 
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We first solve the equations of motion for the surface via power series expansions. Integer 
power series, however, cannot represent nonanalytic functions. Since extremal surfaces in pure 
AdS are nonanalytic at the boundary and the Kasner-AdS soliton is asymptotically locally AdS, we 
expect extremal surfaces in (2.5) to have the same nonanalytic asymptotic behavior. We therefore 
investigate the asymptotic behavior of T(x) and Z(x) to find the corresponding fractional power 
series. 



To probe the asymptotic behavior of T(x), consider the AdS^+i soliton metric (2.4) in Milne 
coordinates: 



/ dz 2 \ 

( (l - z d ) d6 2 + - — — - dt 2 + t 2 dx 2 + dy\ + dy\ 



Implementing a simple coordinate redefinition, we obtain the manifestly static soliton: 

X = t sinh x 



ds 2 



(1 



t cosh x 

v2 



dz 2 \ 
) d6 2 + ^— - - dr 2 + d X 2 + dy\ + dy\\ 



(2.9) 

(2.10) 
(2.11) 



The boundary condition T(x = ±Xq) = i& in the Milne soliton corresponds to a boundary condition 
t{x = ±Xo) = tf, cosh^o in the soliton. The soliton, however, is static, so the entire extremal surface 
rests on the time slice r = tf,coshxo. Using (2.10), we obtain an equation for T(x): 

^Miine^) = tfc cosh x sech x (2-12) 

T(x) is therefore analytic in the Milne soliton. We expect that, for the Kasner-AdS soliton, T(x) 
is analytic as well, and can be expanded: 



T(x) = J2^{x 2 



X J 



(2.13) 



i=0 



where the x 2 term results from T(x) = T(—x) by reflection symmetry of the entangling region. 
The exact solution for x(Z) is known in pure AdS 



±x(Z) 



dzi 



1 



d 



3d -2 Z 2 ^" 1 ) 



2'2(d-l)'2(d-l)' £W 





3d- 2 
2(d-l) 


dr 




2d— 1 




2(d-l) j 





(2.14) 



where 



2xr 



(d-i)r 


2d-l 


2(d-l) 




d 




2(d-l) 





(2.15) 



is the turning point of the surface and 2-F1 is a hypergeometric function. We expand the solution 
in a Taylor series about the point xq and then take the inverse series to obtain an expansion for 
Z(x): 



Z(x) 



(2.16) 



8=0 
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where the Zi are constants. We expect a similar asymptotic behavior in the Kasner-AdS soliton. 



Indeed, substituting the expansions (2.13) and (2.16) into the equations for the extremal surface 



results in a two parameter family of solutions. The remaining two parameters are fixed by requiring 
regularity at the origin. This is difficult to impose starting with local power series expansions about 
±Xo, so we turn to numerical methods. 

One can make the expansion for Z[x) analytic by introducing a new coordinate x defined by 
(xq — x 2 ) d = (1 — x 2 ). The equations of motion remain analytic in x and we have 

oo 

Z(x) = (1 - x 2 ) Y^Zi{l- x 2 Y = (1 - x 2 ) P (x) (2.17) 

i=0 

The factor of (1 — x 2 ) on the right is convenient since then P(±l) is nonzero. Since P(x) is analytic 
everywhere, we can solve for it and T(x) numerically using pseudospectral methods. For a single 
connected extremal surface, this is most conveniently done by taking advantage of the reflection 
symmetry about x = and just solving the equations for x in [0,Xo] with boundary conditions 
T'(0) = 0, Z'(0) = 0,T(xo) = tb and Z(xo) = 0. We have checked our numerics by comparing the 
numerical solutions with known exact solutions in special cases. Our numerical solutions for Z[x) 
in pure AdS and T(x) in the Milne soliton agree with the known solutions to order 10~ 14 . We 
solve the equations of motion for a range of the parameters pi and tj,. By comparing the power 
series solutions with the numerical results, we fix the two free parameters in those solutions. 

The reason that we choose to work with analytical power series rather than the numerical 
solutions is that numerical integration of the surfaces is difficult to perform accurately due to the 
steepness of the surfaces at boundary. In particular, extracting the dependence of the area on 
the cutoff is delicate to do numerically, and much more accurate results can be obtained from the 
power series as follows. We compute the area by substituting the (truncated) series expansions for 



T(x) and Z(x) into (2.7), expand the functional in polynomials, and integrate analytically from 
x = to x = xuv-i where xuv is a symbolic parameter without an assigned numerical value. The 
analytical form of the integral A(x) is a power series in xuv'- 

A(x uv ) = J2 A i(xo-x uv f d (2.18) 

i 

where the Ai are numerical coefficients. In order to determine the area's UV cutoff-dependent 
and independent behavior, the area A(xjjv) must be expressed as a function of Z. To accomplish 
this, we obtain a series expansion for x(Z) by inverting the series Z[x) via a variable change 
(xq — x 2 ) l ^ d = Xq — y 2 , so that we have, as before, an analytic series expansion for Z: 

Z{y) = Y,Z l {xl-y 2 ) 1 (2.19) 

i 

Since this series is analytic, we can easily obtain the inverse series y{Z) and the area A (x (y {Z))) = 
A (Z) as a function of the bulk coordinate Z. A Taylor expansion of A(Z) around Z = yields the 



expected behavior ( 1.3 ). This allows us to immediately read off the leading divergent, subdivergent, 



7 



and constant pieces for a choice of pair {pi, tb). This process is repeated for a representative range 
of pi and tb. We use a least-squares best-fit algorithm to obtain the coefficients of the divergent, 
subdivergent, and UV cutoff- independent pieces as functions of tb and p\. 



3 Results 

3.1 Confinement /deconfinement transition 

We first describe the behavior of the extremal surfaces in the bulkj^] The possibility that extremal 
surfaces will bend toward the singularity and reach it at some nonzero value of the boundary time 
tb may appear to be a concern, as this would make it impossible to evaluate the entanglement 
entropy at any earlier time. We find that this does not occur. The bulk time t increases along the 
extremal surface, reaching a maximum t* at the midpoint x — 0. For tb of order one, the increase 
is small (of order a few percent), but closer to the singularity it becomes larger. For example, for 
Pi = —0.25 and t b = 0.4, t* is 20% larger than tb- 

Next we consider the radial extent of the extremal surface. For a static strip in pure AdSg, the 
extremal surface reaches a maximum value of z given by \27\ 



4r ' 9 / 8 ' C « 1.48£ (3.1) 



vrr(5/8) 



where C is the width of the strip. We find that this formula works remarkably well for the Kasner- 
AdS soliton with C interpreted as the proper width of the entangling region, until z* gets close to 
one. The fact that the cap in the radial direction due to the soliton does not affect the surface 
much when z* < 1 can be understood from the fact that we work in six dimensions, so the effect of 



the cap comes in at 0(z 5 ) in (2.5). It is more surprising that the time dependence does not seem 



to affect the radial extent of the extremal surface (at least in the coordinates used in (2.5)). Eq. 



(3.1) holds even when the extremal surface does not stay near a constant t surface in the bulk. 

When z* does approach one, there is a qualitative change in the extremal surface. Rather than 
one surface connecting the two ends of the strip, the extremal surface splits into two pieces, with 
each one wrapping the cap at z — 1 (see Fig. [TJ. The topology of the extremal surface in the 
x, z, 9 directions changes from K x S l to two copies of IR 2 . Once the extremal surface splits, the 
entanglement entropy becomes independent of the width L. 

This transition is the holographic description of a confinement / deconfinement phase transition 
(HHH]. Entanglement entropy of a confining gauge theory measures the effective degrees of freedom 
between the entangling region and its complement at the energy scale E ~ 4v. The corresponding 
modes of the gauge theory are in the deconfined phase when E is large and the entanglement 
entropy is a function of C. As C grows, one probes progressively lower energy modes. Eventually 



3 Thc specific calculations described in this section were done with xq = 0.12. This is a convenient value since it 
is large enough that the extremal surfaces do not stay just in the asymptotic AdS region, and small enough that 
the extremal surfaces stay connected for a range of tb- It is worth noting that while x$ = 0.12 is a convenient choice 
for the numerical analysis, the qualitative behavior is independent of this choice. 
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z 
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Figure 1: For small proper width C in the Kasner-AdS soliton background, the extremal surfaces remain 
close to the boundary and resemble surfaces in pure AdS. The area of the surface depends on C, so the 
entanglement entropy is a function of C and the modes which contribute are deconfined. As the length 
of the entangling region becomes larger, extremal surfaces extend further into the bulk. The soliton cap 
begins to affect their geometry, and they start leveling out. For some critical length, the surface splits into 
two separate surfaces (displayed above in dashed blue). At this point, the area of the extremal surface is 
independent of C, effectively signaling that 1/C is below the mass gap of the confined gauge theory. 

one hits the mass gap for the confining theory and there are no modes of lower energy. When this 
happens, the entanglement entropy becomes independent of C. 

The sign of pi clearly plays an important role in determining when this transition occurs. 
When pi < 0, the proper length C grows as one approaches the singularity and the confine- 
ment/deconfmement transition takes place at early time. When p\ > 0, C grows as one moves 
away from the singularity, so the confinement/deconflnement transition occurs at late time (see 



3.2 Entanglement entropy 

We now proceed to present the behavior of the area of the extremal surface as a function of tj, 
and p\. As noted above, this is divergent, so we introduce a small cutoff at Z = e and expand 
A[Z = e). It is convenient to factor out the trivial integral over yx,yz, 9 and write: 



As expected, the subdivergent piece 1/e 2 vanishes, and the leading coefficient is proportional to 
the area of the entangling region: 



Fig.g. 




(3.2) 



1 A 



bdy 



2 



-P2+P3 _ 




1-P1 



(3.3) 



C 3 = o 



3LiL 
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(a) (b) 

Figure 2: (a): The maximum value of z as a function of tb, the boundary time. Each line corresponds 
to a different value of p\ G [—0.3,1]: surfaces with positive p\ shrink towards the boundary as tb — > 
while surfaces with p\ < approach the cap. Surfaces with p\ = maintain constant z*. (b): A plot of 
z*(tb) = 1-48/3 for small tft and p\ 6 [— g,l]. This describes the radial extent of minimal surfaces until 
~ 1 when the surface splits into two, signaling a confinement /deconfinement transition. 



where the factor of 2 is due to the fact that the yi and 9 are maximally extended, so &1Z consists 
of two identical disconnected pieces. The last equality comes from the Kasner conditions (2.2). 
The factor of 1/3 in the first equality can be understood geometrically: the bulk is asymptotically 
locally AdSg, so near Z = 0, the area functional is approximately: 

A AdS = j ^dZ (3.4) 

The Z 4 factor results in a 1/3 contribution to the integral. This result moreover agrees with [5], 
where the computation is done for the pure soliton for an entangling region identical to 1Z (tb, Xq). 
We next turn to the subleading divergence. We find Ci(i&) is a simple power law: 

Ci(tb)=l(Pi)C l} 

7 (Pi) = a + bpi + c Pi ( 3 - 5 ) 
/3 (pi) = -1 - pi 

where a, b, and c are numerical constants. 

Considerations of covariance and dimensional analysis lead us to conclude that the C\ coefficient 
is an integral of a dimension two geometric quantity. This can be the scalar curvature of the 
background spacetime, the induced metric on the boundary of the entangling region, or a scalar 
constructed from the square of the extrinsic curvature, K^ u p . Both scalar curvatures vanish in our 
case, so the two possible contributions are: 



A hdy ( 7l i^ p Kl o + l2 K, u 'KV P ) = -^ L [li(l - Pi) 2 + 7 2 (1 - pi)} (3.6) 



for some constants 71 and 72 determined by the field theory. We therefore expect: 

Cl = 2t- x - pi [-( 7l + 72 ) + 2 Pl7l + ( 72 - 7l )p?] (3.7) 
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The power of t& in the time dependence of this expected behavior matches our results, (3(p 



-1 — pi. We may calculate 71 and 72 by setting the quadratic prefactors in (3.5) and (3.7) equal: 

1 



(a + bpi + cpf) = -(71 + 72) + 2pi7i + (72 - 71)^1 (3. 



6 

7i 



V &\ If ^ 



In other words: 



72 2 V" 2 V" ' 2 

This is an overdetermined system and has no solution unless a = — (b + c). We find that our 
coefficients indeed satisfy this relation and furthermore determine 71 « —0.1032 and 72 ~ 0.0238. 
To the best of our numerical accuracy, these constants have a simple ratio — = — y. This is 
sufficient to determine the extrinsic curvature contributions to the 1/e divergence in general for a 
six dimensional bulk. 

We now discuss the UV finite part of the entanglement entropy, Co(tb). This is always negative, 
but has several qualitatively different behaviors depending on the value of p\. Recall that for a 
static strip in pure AdS6, cq behaves as |24| : 

Co, static 3^" (3.10) 

A naive generalization of this for a time-dependent metric would be: 

cooc^f (3.11) 
where, as before, L is the proper width of the entangling region. We expect to find agreement with 



(3.11) for the Milne soliton (pi = 1) since this is static by the coordinate transformation (2.10). 
Indeed, we find that for p\ between 1 and approximately 1/2, co(£&) is well fit by a simple power 
law 

^r ?1 = |r^ (3.12) 

The agreement with this power law is illustrated in Fig. [3j We can attribute this agreement to 
the fact that for 1/2 < p\ < 1, the time dependence in the x direction dominates over the time 
dependence in the yi directions. 



If Co obeyed (3.12) for all pi, there would be an abrupt change in the behavior of Cq at the 
critical value pf* = 1/4: surfaces with p\ > p^ would have cq diverging at early tf, and going to 
zero at late while surfaces with p\ < p± would have Cq vanishing at earlier times and diverging 
at late times. This is not what we find. Instead, for p\ < 0.5, Co(t&) is no longer fit by a simple 
power law. Fig. [4] shows the behavior of cq for small, positive p\. For 1/4 < p\ < 1/2, cq is less 



steep at early times than the static case would predict, while still diverging. For < p± < 1/4, 



|co(t&)| decreases as one moves toward the singularity as suggested by (3.12), but eventually turns 



around and starts to increase. We believe that it will eventually diverge as tf, — > 0, although it 
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Pi 



Figure 3: A plot of C(pi), where co(t) oc 4 , f° r | < Pi < !■ The best-fit line is C(pi) = 1 _ 4pi. 
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Figure 4: cq as a function of time for < p\ < 1/2. As p\ decreases, the cq curves level out, while still 
apparently diverging at early times, in contrast with the static prediction, which would imply that they 
go to zero for p± < 1/4. 



is hard to numerically find the extremal surfaces very close to the singularity (and one cannot 
trust the classical bulk solution very close to the singularity anyway). For all of these cases, the 
width of the entangling region decreases at early times and the extremal surfaces stay close to 
the boundary. The fact that |co(t&)| diverges at early time likely stems from the fact that it is 
measuring entanglement of modes with shorter and shorter wavelength. 

It is clear from Fig. [i]that all the curves cross at « 0.95. In other words, Co(0.95) is essentially 
independent of p\. This is a result of the fact that for t ~ 1, the Kasner metric is independent of 
Pi, and hence the same is true for the Kasner- AdS soliton. Since the extremal surfaces are staying 
close to a constant t surface when t m 1, their area is essentially independent of pi. 

Finally, for p\ < 0, |cfo(t&)| monotonically decreases at early times as illustrated in Fig. [5} In this 
case, the proper width of the entangling region grows and the extremal surface probes progressively 
deeper into the bulk. This is expected to continue until the confinement / deconfinement transition 
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Figure 5: Co as a function of time for —0.3 < p\ < 0. The curves correspond to pi = —0.03, —0.11, 
-0.17,, -0.22, -0.26, and -0.30. 

occurs, when the extremal surface hits the cap of the soliton and splits into two disconnected 
surfaces. 

To summarize, we have computed the entanglement entropy of a confining gauge theory in 
a background spacetime which is a product of Kasner and a circle. In addition to the leading 
divergence which is proportional to the area, we found a subleading divergence proportional to the 
square of the extrinsic curvature of the entangling region. Near the singularity, we find that the 
behavior of the UV finite part of the entropy, cq, can differ significantly from the naive extrap- 
olation of the static result: the behavior of the width of our entangling strip is more important 
than the behavior of the area of its boundary. If the width goes to zero at early time, then cq 
diverges, while if it grows, then |co(t&)| decreases near the singularity and eventually reaches a 
connnement/deconfinement transition. This shows that the approach to the singularity can be 
used to determine the energy scale at which a confinement / deconfinement phase transition occurs 
in a gapped field theory. 

In cosmology, one often considers states with nonzero temperature, and a confinement/ decon- 
finement transition occurs at a time when the temperature reaches a critical value. In contrast, 
we are not working with a state with nonzero temperature. Our state is given implicitly by our 
choice of bulk solution (since different states in the dual field theory correspond to different bulk 
geometries). However, one cannot describe this state as a vacuum state since the time dependence 
of the background causes particle creation. A vacuum state at one moment of time will not remain 
the vacuum at later time. This particle creation is contributing to Co(i&). 

It would be interesting to compare our calculation of the entanglement entropy, which is ap- 
propriate for a strongly coupled confining theory, with the entanglement entropy of a free field in 
Kasner (see, e.g., p8j). 
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